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I. INTRODUCTION 



Recent work on the dynamics of wave-particle interaction has led to extensive use of N- 
body hamiltonian models in parallel with the more traditional kinetic approach. The present 
paper aims at discussing to what extent the two approaches agree in the limit iV — > oo, 
where iV-body dynamics formally reduces to kinetic theory. This is a classical problem of 
statistical physics, notoriously unsolved for particles interacting through short-range forces, 
where it amounts to deriving the Boltzmann equation from the Liouville equation : sys- 
tematic rigorous derivations of the kinetic equation from BBGKY hierarchy are still lacking 
- notwithstanding the pioneering work of Lanford and King, limited to short timescales 
1^ , 20] , and recent advances However, for long-range forces, and more precisely for 



smooth enough mean-field interactions, the formal limit N — > oo commutes with the dynam 
ics 



2^j2q| . We show in this paper how the mean-field methods apply also to wave-particle 
interactions. 

A physical motivation for this work is that wave-particle interacting systems are typical 
of plasmas and common to many physical phenomena. The paradigm of such interactions is 
provided by the self-consistent hamiltonian H^.' describing the evolution of iV particles and 
M Langmuir waves [p^j5|,P- p^fi^ , [^j23] , |2^j27| , |28[] . In particular this hamiltonian enables a 
mechanical approach of classical plasma problems like Landau damping and beam-plasma 
instability, by treating Langmuir waves as M harmonic oscillators self-consistently coupled 
to iV quasiresonant beam particles. For simplicity we present our results in one space 
dimension with periodic boundary conditions, which conforms to the physical conditions 



considered in models of plasmas fl J12|J15 l. 

The basic characteristic of such models is that particles do not interact directly with 
each other : they only interact with the modes ; symmetrically, the modes do not interact 
directly with each other : they only interact directly with the particles. Inasmuch the modes 
are spatial Fourier components of some fields, these components are not localized spatially : 
this invites to describe the many-body limit iV — > oo as a mean-field limit and enables us 
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to apply the techniques which succeed in the case of particle-particle mean-field coupling. 

The present work takes advantage of this observation to show that the kinetic limit 
iV —>■ oo and the time evolution over any time interval [0, T] commute. Our result implies 
that numerical simulations with increasing number of particles behave ever closer to the 
predictions of kinetic theory (if one uses N 'large enough'...). 

A preliminary form of our result was annouced in reference [ 14| . In Sec. II we describe 



the model and its evolution equations. The main results are stated in Sec. III. Sec. IV is 
devoted mainly to a finite N estimate and a technical remark, preparing the proof presented 
in Sec. V. The final section is devoted to the conclusion. 

II. SELF-CONSISTENT HAMILTONIAN AND KINETIC LIMIT 

We consider a system of TV particles with respectively mass m r , charge q r , position x r 
and momentum p r , interacting with M waves with respectively natural frequency u)j , phase 
8j and intensity Ij. The evolution of this system is described by the hamiltonian 

N 2 M N M 

H sc M = E -tr + £ ^ - £ EE qrkffrfiTj C os(k jXr - e 3 ) (i) 

r=l Zm r j=l r =l j=l 

where the first term corresponds to free particles, the second term to free waves (harmonic 
oscillators) and the third term to their coupling. The coupling constants are expressed in 
such a way to ease the kinetic limit N —>■ oo : we shall keep the 'wave susceptibilities' (3j 
constant in this limit. A simple change of variables enables one to ensure that all coefficients 
(3j > 0, which is assumed in the following. The overall coupling factor e in the interaction 
term of (P emphasizes our interest in the weak-coupling regime (e 1) P,f["lJ|. 

Assuming periodic boundary conditions, the particles move on (M/L) = Sl and the 
wavenumbers are quantized (kj = rij27r/L for some integer nj). The phase space of this 
system is thus (Sl x M) n x Z m where Z = x M + for each mode. 

The natural scaling of our model in the limit A — * oo is easily deduced from its equilib- 
rium (Gibbs) thermodynamics |TB|. Then the energy E = H sc and the wave intensities Ij are 
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extensive (i.e. 0(N)), and the coupling constant scales as e = OiN- 1 / 2 ). The extensivity of 
wave intensities can easily be interpreted as, in the physical regime of the model, we expect 
particles to be mostly resonant with the waves, each such particle contributing then to wave 
intensities by evolving in their potential well. This prompts us to introduce intensive wave 
variables 



z 



N-^Zj = N" 1 ' 2 JWje-** = \zj\e~^ (2) 



for which Z reduces to C, and renormalized coupling constants /3'j = N 1 / 2 e/3j. 
The evolution equations of the hamiltonian ([I]) read 



x r = p r /m r (3) 

• M 

Pr = - z*e~ ik ^) (4) 

A 3=1 

i N 

Zj = -tuj j0 z, + -pfa 1 £ Ire tk ^ (5) 



r=l 

To simplify calculations, we introduce non-canonical variables, namely particle velocities 



v r = p r /m r (6) 



and mode envelopes 



a, = z 3 e^ (7) 



bringing to the form 



x r = v r (8) 
ia M 

y r = _1!L J2 p' j (a j j k i x '~ iu i ot - a* j e- ik * Xr+iu * ot ) (9) 

r j=\ 

i N 

a- = —p-kj 1 q r e- ik ^ r+l ^ ot (10) 

The usual space of kinetic theory is Boltzmann's /i-space A = Sl x M. The positions 
and velocities of the N particles determine a distribution / on A : 
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1 N 

/(*, = Jj E - *r(*) - »V(t)) (11) 
r=l 

which is normalized to unity (/ A f(x, v, t)dxdv = 1) irrespective of the number N of particles 
(for simplicity, we assume a single species : all q r = q > 0, m r = m). The kinetic limit, 
formally iV — > oo, corresponds to considering a sequence of A^-particle distributions / 
converging to a distribution f°° in the weak sense for a natural space of observables T>. 
Denote by T the space of positive normalized distributions on A with finite momentum and 
kinetic energy, i.e. T = {/ G -^ 1 (A, dxdv) : / > 0, / fdxdv — 1, / v 2 fdxdv < oo} and define 
on T the bounded-Lipshitz distance 

dbL (/,/') = SU P I / 4> fdxdv - I <j)f'dxdv\ (12) 

with the set of bounded, Lipschitz-continuous normalized observables 

V = {0 : A - [0, 1], |0(x, v) - 0(x', t/) I < || (x, v) - (x', v') || V(x, u), (x', «') G A} (13) 

Here A is equipped with the distance ||(x, v) — (x',v')\\ = a(\(x — x')modL\ + r\v — v'\), 
where a -1 and r are respectively convenient length and time scales to be chosen below. 
Then we consider the distance on Z M , 

M 

\\a — a'\\ — y^Wjlcij — a'-| (14) 

3=1 



where real positive coefficients Wj will be chosen below in ([31]) , and \aj\ is the modulus of 
the complex number a,-. Our distance on JF x Z M is just 

||(/,o)-(/' ) a / )|| = d Ui (/ ) / / ) + l|a-a / |l (15) 

The kinetic evolution system of equations, dual to (p|-p|-p!0D, is the system 

ia M 

d t f + vd x f + —J2( 3 j (a 3 e ik ^- l ^ ot - a*e-^ x+i ^ ot )d v f = (16) 
2m i=1 



dj = iqp-kj 1 f(x, v, t)e-^ x+ ^ ot dxdv (17) 
This dynamics leaves T x Z M invariant. 



III. MAIN RESULTS 



The self-consistent dynamics preserves two constants of the motion, namely total 

energy H and total momentum P = J2 r Pr + J2j kjlj- 111 the kinetic limit, we consider the 
normalized constants h = H/N and p = P/N : 

2 I 1 2 

KM = I ^-fdxdv + ^oH^- [ ^qkj^'^ia^-^fdxdv (18) 

JAZ j Z J A j 

I 1 2 

p(/, a) = / mvfdxdv + J2 k j^T~ ( 19 ) 
Ja j 2 

where 3? denotes the real part. For any finite iV and h, the energy surface H^ C ,M = Nh in 
x Z M is compact, and the vector field is continuous and bounded on it. This 

ensures that the dynamics generates a group for all initial conditions. 

Moreover, the first variation of the dynamics (]8]-|9|-|iT)|) generates a linear operator M. = 
d(x r ,v r ,dj)/d(x r ,v r ,aj), depending continuously on (x r ,v r ,a,j). As the energy surface is 
compact for any given N, M. is bounded. With the specific form of H sc , we show that, with 
appropriate choice of the constants Wj : 

\\M || < r- 1 + 7 [a(t)] (20) 

where 

T =(-£#y 1/s . (2D 

7[aW] = E^ilaiWI ( 22 ) 
j 

The positive function 7[a(s)] is continuous on the energy surface, on which it has an upper 

bound uniform with respect to N. 

The kinetic limit, — ► oo, admits a similar bound, ensuring the existence and uniqueness 
Theorem : Given initial data (/ , a(0)), (/q, a'(0)) G T x Z M , with h = h(f , a(0)) and 

h' = h(fb, a'(0)), the kinetic evolution equations ( ]TB| - |r7D generate for all times t > states 

(ft,a(t)) and (fl,a'(t)) respectively from these data. Moreover, 
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d bL (ft, ft) + Ht) - a'(t)\\ < e ct (d bL (f , f ) + \\a(0) - a'(0)||) (23) 

for some C = C(h , h' ) < oo. 
This theorem implies the 

Corollary : Given a distribution /g° G T and a sequence of finite- iV Dirac distributions 
f^^T for particle initial data, such that lirn/v-+oo ^ftiX/o^ /q°) = 0, given initial waves 
a(0) G Z M , and given any time T > 0, consider for all < t < T the resulting distribu- 
tions //^ and waves a N (t) generated by H^ C ,M and the kinetic solution (/ t °°, a°°(t)). Then 
limjy^oo dbL^ft 1 , ft°) = and lim^^oo a N (t) = a°°(t), uniformly on [0; T]. 
In other words, the following diagram commutes for all t > : 
(f N ,a(0)) (8 "1P (f t N ,a N (t)) 

I iV -> oo | iV -> oo (24) 



IV. PRELIMINARY REMARKS 

For given iV and finite energy H, the first variation M. of the dynamics (@)-(@)~(H3) nas 
bounded norm (with the L x distance) : 

N M 

\\(5x,Sv,5a)\\i = N" 1 a(|&r r | + r\Sv r \) + ^w,-|<%| (25) 

r=l j=l 

N M 

= N- 1 Y,\\(8i r ,6v r )\\ + Y,v>i\foj\ 

r=l j=l 

We readily find 

\6x r \ = \5v r \, (26) 

N N 

\5aj\ = N- l \ ^2 P' jq e ik ^-^ ot 6x r \ < iV" 1 /^^ \Sx r \ (27) 

and 



S 



M qft 



\8v r \ = |£ m&(e ,k * x '-' u * ,t (a j k j 6x r - i8aj))\ (28) 

3=1 m 

M q&- M q&- 

<E^%K-I-M + E^-KI ( 2 9) 

so that 



j=1 m p{ m 



N 



\\(8x,6v,8a)\\i < N-^a^ir^ait)} ■ \5x r \ + r\5v r \) 

r=l 

M N M Q QI 

+ J2 w i N ~ 1 f 3 'j ( lYl \5x r \+aTj2 — - \Saj\ (30) 

3=1 ' r=l j=l 171 

with 7[a(s)] defined by (|22|) . The four causes for the divergence of trajectories in x Z M 
are saddle points (in (x,v) plane) associated with maxima of the modes' potentials (the \aj\ 
contribution to 7[a(t)]), velocity shear (the velocity term), the dependence of the modes 
source on the particle positions, and the dependence of the saddle points themselves on the 
mode envelopes. 

An appropriate choice of constants a, r, Wj keeps the estimates as small as possible. 
Thus let 

Wj = awoPp (31) 

and solve 

r-^^oE^f^— (32) 
i mw o 



This leads to the expression of r announced in fl2T ) and to 



w,= {qm)-UY^(3fY\ (33) 



3 
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so that (|30|) reduces to 

|| (8x, 6v, 8a) || < r- 1 1| (6x, 8v, 8a) || + 7 [a(t)]^ 1 ^ a\8x r \ (34) 



which implies (p0|). Constant a remains arbitrary, as it only determines the scale of the 
distances in A and Z, and (PD]) is homogeneous (degree 1). Considering only the restricted 
dynamics on A, with 8a = 0, (|3~0"D straightforwardly leads to the continuity equation 



\\(Sx,5v)\\ <^[a(t)]\\(5x,5v)\\ (35) 

with 

7 7 [a(t)] = max(r" 1 ,7[a(t)]). (36) 

Note that 7'[a(t)] is bounded uniformly in time, as the positive function 7 [a] is bounded 
above on the energy surface by a function which does not grow faster than h 1 ^ 2 in the large 
energy limit. More precisely, let A > solve A 2 J2j ^Jo (iP'jkj I m ) 2 = 2/i + J^j^joiqPj/kj) 2 . 
Then J2j \0jkj(ijq/m\ < (q 2 /m) J2j ^joft'/i^- + k 2 \/m). It should be noted once more that 
(|35l) reflects that the divergence rate in (x, v) is controlled by velocity shear and by saddle 
points of the pendulum-like potential depending on wave amplitudes. The latter situation 
typically corresponds to a trapping regime for large enough wave intensities. 
Finally, note the following 

Proposition 1 : Let Y : A — > A be a Lipschitz mapping with constant L > 1 on A, and 
/i, v G T . Then : 

sup I / o Yd(fi — v)\ < Ld bL (fi, v) (37) 
4>ev JA 

Proof : Clearly L _1 o Y E V for any <\> E V. Hence sup^ eX > I Ia^ 1( P Yd(ji — u)\ < 



V. PROOF OF THE MAIN RESULT 

The proof of theorem 1 uses the fact that the two types of degrees of freedom have no 
'self-interaction. Indeed the motion of particle r is completely determined by its initial 
position and velocity and by the modes history, i.e. the data of the modes %(.) over a time 
interval [s, t] defines the vector field G so that: 

j(x r {t),v r {t)) = G[a(t)](x r (t),v r (t)) (38) 



This vector field is Lipschitz-continuous on A according to fl3"5]) and subsequent remarks. 



Thus Cauchy-Lipschitz theorem ensures the existence and unicity of the flow T: 
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(x r (t),v r (t)) = T t!S [a(.)](x r (s),v r (s)) (39) 
By duality the measure fi s on A is transported by the flow to 

Ht = H.oT tt t[a(-)] (40) 

Similarly, the evolution of mode j is also completely determined by its initial data aj(s) and 
by the history of the measure on particle phase space A in the right hand side of ( |l~7|) which 
defines a flow S by 

a,j(t) = S t , s \jJL.]aj(s) (41) 

Solving kinetic equations ([16| ) -([T7| ) with initial data (//q, a(0)) amounts to finding a fixed 
point of the coupled system ((40]) - ( (41]) in the space {T x Z M ) R . Our strategy now follows 



that of Neunzert p2| and Spohn [26fl, who considered direct particle-particle interaction of 
mean-field type. 

Thus consider two solutions (/., a(.)) and (/', a'(.)) of fllED-fll?!). To shorten notations, 
write b = a' and denote by d\i = fdxdv and du = f'dxdv the corresponding measures. Their 
distance at time t satisfies 

\\(jit,a{t)) - (vt,b(t))\\ = 4(/iooT 0)t [ fl (.)], l /ooT Q , t [6(.)]) + ||^,o[//.]a(0) - St,oM&(0)|| (42) 
where 

||£ t ,oHa(0)-S t >.]&(0)|| <d 1 (t) + d 2 (t) (43) 



4l(a*o ° 7o,t [«(•)]> ° T 0) t[b(.)]) < d 3 (t) + diit) 



(44) 



and 



d 1 (t) = \\S t , Q {>M]a{0)-S t , o lpL.)b(0)\\ 



(45) 



4(t) = ||^o[ / u.]KO)-^,oH6(0)|| 
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(46) 



d 3 {t) = d bL (n o T 0>t [a(.)], is o T 0>t [a{.)]) 



d A (t) = d bL (v o T 0>t [a(.)], vq o T 0>t [b(.)]) 

Straightforward integration of ( |TTD shows that 

d^t) =d 1 (0) = ||o(0)- 6(0)|| 

because the flow S[fi] is just a translation in Z M . 

To estimate d^ we integrate (|T7|) with the right hand sides given by fi, and v. 

d 2 (t) = Y,w j qP'k- 1 / / e- lk > x+i ^ s d((i s - v s )ds 
■ =1 Jo Ja 



< \PIt~ x I dbd^s, v s )ds 
Jo 



d(n s - u s )ds 



In fl52|) the inequality uses the fact that a(l+cos(kjX—8)) / kj G T> and a(l+sm(kjX- 
T> for any real 8, provided that 2a < min kj. 

Estimating d^ is also straightforward, as proposition 1 implies 

d 3 (t) < d 31 (t)d bL (fi , v ) 

provided that d^i{t) is a Lipschitz constant for T( i o[a(.)]. Now, V(x,t>), (x',v') G A, 

\\T t>0 [a(.)](x,v)-T ti0 [a(-W^')\\ 
< \\(x,v) - (x',v')\\ + J* \\G[a(s)]T S)0 [a}(x,v) - G[a(s)]T s , [a](x' , v')\\ds 

< \\(x,v) - (x',v')\\ +f* 1 '[a(s)}\\T S!0 [a}(x,v) - T s>0 [a](x' ,v')\\ds 
Hence, dsi(t) < 1 + J * / y'[a(s)]d 3 i(s)ds, which implies 



0^31 (t) < exp / 'y'[a(s)]ds 
Jo 



by Gronwall's lemma. 
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Finally, 



d 4 (t) = sup ((f) o T t ,o[a(.)] -<j>o T tfi [b(.)])du < d m (t) 
</>e£> J A 



where 



d 40 (t) := sup||T^ [a(.)](a;,«) - T t0 [b(.))(x, v)\\ < d 41 (t) + d 42 (t) 

A 



with 



d 4 i(t) := sup || / (G[a(s)]T S:0 [a(.)](x,v) - G[a(s))T s0 [b(.)](x,v))ds\ 
a Jo 

< / j'[a(s)]d 40 (s)ds 



(57) 



(58) 



(59) 
(60) 



d 42 (t) := sup || / (G[a(s)]T sfi [b(.)](x,v) 
A -Jo 



-G[b(s))T sfi [b{.))(x,v))ds\ 



< 



/ sup||G[a(s)] - G[b(s)]\\ds 

JO A 



Definition (13~8T) shows that 



(61) 
(62) 



\\G[a(s)}(x,v)-G[b(s)}(x,v)\\ = |^ £ ^{a 3 {s) - 6,( S ))e^| < r" 1 ||a(a) - 6(a) || (63) 



m 



Now define <p(t), a majorant of the sum d\(t) + ds(t), and ^(t), a majorant of the sum 
d 2 {t) + d 4 (t) as 



<p(t) = \\a(0) - 6(0)|| + e^'' a (* s 4(ft, v ) 



(64) 



d 5 (t) = d 2 (t) + d 41 (t) + d 42 (t) 

Then previous inequalities for the d{ lead to 

d 5 (t) < \plr~ l [ d bL (fi s ,u s )ds+ [ 'j'[a(s)]d 5 (s)ds + T~ 1 [ di(s)d& 
Jo Jo Jo 

and dbL,((J> s , u s) + di(s) < d\(s) + d 3 (s) + d 4 (s) < <p(s) + d 5 (s) so that 

d 5 (t)< f V2r- 1 (p(s)ds+ [\V2r- 1 +j'[a(s)])d 5 (s)ds 
Jo Jo 
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(65) 



(66) 



(67) 



which Gronwall's inequality readily estimates by 

d 5 (t)< f ' y/2r- 1 (p(s)e^ T ~ 1+ ^ la(u)])du ds (68) 
J o 

The resulting complete estimate 

\\{fi t ,a(t))-(u t ,b(t))\\<^(t)+d B (t) (69) 

depends on two functions 7'[a(s)] and <p(s). Note that <p(Q) = ||(/xo,a(0)) — (i/o,£>(0))|| and 
^5(0) = 0. Estimate (|69| ) does not grow faster than exponentially, with upper bound on its 
growth rate 

C = V2T- 1 + 2 sup j'[a(s)] (70) 

0<s<t 

which is bounded by a function of h as discussed in Section IV. This completes the proof 
of the theorem. 

The corollary follows in a standard way. 

Remark : our estimate for the growth rate C in the kinetic case is larger than the finite- N 
estimate for \j\4\ in phase space. This is due to the fact that the distance dbi makes no 
distinction between x-components and f-components, while estimates of Sec. IV relied on 



treating these components of the phase space points separately to obtain fl30|) . 



VI. CONCLUSION 



This work supports theoretically the use of full iV-body dynamical schemes [|||| |lT] , |r6 -18 



to study the wave-particle interactions, as an alternative to kinetic-theory based models. 
However the regularity of the limit N — ► 00 is tempered by the rapid growth of the right 
hand side in the upper bound (^3|) . 

It also identifies the fundamental cause of phase space mixing and approach to equi- 
librium in this many-body system : particles passing near the instantaneous saddle points 
associated with the modes undergo exponential dichotomy, with a divergence rate controlled 
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by amplitudes \zj\ = |aJ. This implies that the phase space regions where discrepancies be- 
tween the kinetic description and the finite- N description show up most rapidly correspond 
to the neighbourhood of the 'separatrices' associated with the envelopes in the particles' 
//-space A, as was observed in numerical simulations for M = 1 by Guyomarc'h P7|JI^| . 
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